Magnetic susceptibility due to disorder-induced 
neutral solitons in interacting polymer chains 



Marc Thilo Figge* Maxim Mostovoy, and Jasper Knoester 
Institute for Theoretical Physics and Materials Science Center 
University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands 

February 1, 2008 



Abstract 

We study the magnetic response due to neutral solitons induced by dis- 
order in polymer materials. We account for interchain interactions, which, 
if sufficiently strong, result in a bond-ordered phase, in which the neutral 
solitons are bound into pairs. We analytically calculate the corresponding 
pair size distribution. As the spins of the solitons have a distance dependent 
antiferromagnetic coupling, this allows us to calculate the magnetic suscep- 
tibility in the ordered phase. At low temperatures, the result deviates from 
the usual Curie behavior in a way that depends on the relative strength of 
the disorder and the interchain interactions. We compare our results to the 
observed magnetic susceptibility of irans-polyacetylene and we suggest new 
experiments extending towards lower temperatures. 
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1 Introduction 



In recent papers we have studied the effect of off-diagonal disorder on the lattice con- 
figuration of a half-filled Peierls-Hubbard chain with a doubly degenerate ground state 
H|, [|. An example is the 7r-conjugated polymer trans-polyacetylene, which in the absence 
of disorder has a uniformly dimerized ground state. As the energy does not depend on 
the sign of the lattice dimerization (which determines whether even or odd bonds are 
short), this ground state is doubly degenerate. For a single Peierls-Hubbard chain we 
showed that arbitrarily weak off-diagonal disorder induces neutral solitons (kinks) in the 
lattice dimerization, which interpolate between the two degenerate bond alternations. 
Even though the creation energy of a neutral soliton is rather large (of the order of the 
gap in irans-polyacetylene), the energy loss is compensated by allowing the sign of the 
dimerization to adjust to the electronic disorder fluctuations. Off-diagonal disorder in 
conjugated polymers (i.e., disorder in the hopping amplitudes of the 7r-electrons) may 
originate from random chain twists, which decrease the overlap between the 7r-orbitals of 
neighboring carbon atoms. The density of disorder-induced neutral solitons in a single 
chain is proportional to the disorder strength [[[], ||. 

In £rans-polyacetylene the disorder strength is presumably rather large, as the average 
conjugation length in this material is of the order of several tens of carbon atoms only||. 
One would then expect the neutral solitons to contribute significantly to the polymer's 
magnetic and optical properties, as they carry spin | and result in the appearance of 
electronic states inside the Peierls gap. 

There is, however, no direct evidence for the existence of a high density of solitons 
in undoped trans-polyacetylene. Electron spin resonance (ESR) experiments report only 
about one free spin per 3000 carbon atoms ||, |J |, ^]. Moreover, it appears difficult 
to observe neutral solitons in optical absorption experiments, as, contrary to what is 
expected from the Su-Schrieffer-Heeger (SSH) model ||, they seem not to give rise to a 
clear midgap absorption peak ||. This may be explained by assuming that the on-site 
Coulomb repulsion U is strong enough to shift the midgap peak towards the absorption 
edge resulting from interband transitions flO], [H| . Since both the peak and the absorption 
edge are significantly broadened by the quantum lattice motion |T2| , there may be no 
clear distinction between them. 

X-ray scattering data do yield some indirect evidence for disorder-induced kinks. There 
still is considerable disagreement whether neighboring carbon chains are dimerized in 
phase (P2 1 /a space group) or in anti-phase (P2 1 /n space group) O, ITJJ. While this 



may originate from different preparation methods leading to different space groups, it 
has also been pointed out that the disagreement may result from a high density (of the 
order of several percent) of kinks that locally change the relative sign of the dimerization 



in neighboring chains pL6| . Yet, it is not clear how such random changes would lead 
to sharp peaks in the x-ray spectra. 

From the above it appears that a clear signature of the effect of disorder-induced 
solitons is still to be found. This has motivated us to study the magnetic response of 
disorder-induced solitons in more detail. A proper modeling of the magnetic susceptibility 
involves more than a calculation of the density of solitons in an isolated chain. As 
we noted in Refs. |], 0, the actual density of neutral solitons (and thus of spins) is 
determined by the competition between disorder and interchain interactions, as the latter 
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lead to confinement of soliton-antisoliton pairs and may restore the long-range bond 
order. Moreover, at low temperature, the exchange between the spins of neighboring 
solitons on a single chain tends to bind them into a singlet state WJ§, which has no 
magnetic response. Thus, the magnetic susceptibility and, in particular, its temperature 
dependence should be expected to depend strongly on the interplay between disorder and 
interchain interactions. 

In this paper, we focus on the magnetic susceptibility of disorder-induced solitons in 
the phase with long-range bond order. In this case, solitons occur in isolated pairs of 
random size dictated by the disorder realization. We study the statistics of these pairs 
by mapping the problem on the anisotropic random-field Ising model, which is treated 
in the chain-mean-field approximation [ft"8|| . This mapping is analogous to what we did 
in Refs. [|l], |2| for isolated chains. We briefly explain the mapping in Sec. |2] and discuss 
the phase diagram of this model. In Sec. we express the magnetic susceptibility of 
an ensemble of soliton-antisoliton pairs in terms of the, as yet unknown, distribution 
of exchange interactions. In Sec. [|, we calculate the distribution of soliton-pair sizes 
using the saddle-point method. From the pair size distribution we derive in Sec. |5] the 
distribution of exchange constants. The latter is used in Sec. ^| to calculate the magnetic 
susceptibility of trans-polyacetylene. We find that the low-temperature behavior of this 
susceptibility deviates from the Curie law and we fit our results to the experimental data 
obtained in Ref. [[/J. In Sec. [7|, we summarize and conclude. 



2 Solitons in Interacting Disordered Peierls Chains 

In Ref. [0 , we have shown that the statistics of neutral solitons in isolated weakly disor- 
dered Peierls chains can be studied using the one-dimensional random-field Ising model 
(RFIM). In this mapping, Ising variables a m = ±1 (m = 1, M) are defined on the sites 
of a lattice with lattice constant d. These variables play the role of the sign of the dimer- 
ization in the Peierls chain, while the random "magnetic" field h m at site m represents 
the off-diagonal disorder, which locally lifts the degeneracy between the two dimeriza- 
tion phases in the Peierls chain. Two neighboring sites on the lattice having different 
Ising variable, correspond to the occurrence of a soliton in the Peierls chain. There- 
fore, the creation energy \i of a soliton in the Peierls chain is equivalent to the exchange 
interaction between neighboring Ising spins. In the SSH model of trans-polyacetylene 
H = 2A /7r ~ 0.5 eV (A is the dimerization). We emphasize, however, that this map- 
ping is not limited to the SSH model, but also holds in the presence of electron-electron 



interactions, in which case the value of \x is smaller [jT^, [20|j . 

Our approach may easily be extended to account for three-dimensional effects: in- 
terchain interactions (electron hopping, eleastic forces, or Coulomb interactions) tend 
to favor a coherence of the dimerization pattern on neighboring chains, which in Ising 
language translates into an interaction, 2W, between spins on neighboring chains. As 
for quasi-one- dimensional materials, like conjugated polymers, W <C fi, we are thus deal- 
ing with a strongly anisotropic random- field Ising model [pjfl. The anisotropy allows 



one to treat the interchain interactions in a mean-field way, an approach known as the 
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chain-mean-field approximation. The energy of the resulting Ising model is given by: 



E{a m } = £ 



M V 



m=l 



7j/(l ~~ Vm+i) ~~ h m a m — Ba r< 



Here, the first term describes the energy cost for creating |l]m=i(l — ^m^m+i) kinks 
and the second term describes the interaction energy with the random magnetic field. 
The latter is assumed to have a Gaussian distribution with zero mean ((h m ) = 0) and 
correlator: 

(hmfon) £ &m,n • (2) 

where e is the disorder strength. We consider the case of weak disorder: e -C /i 2 . Finally, 
the third term in Eq. ([]]) describes the interchain interactions in the mean-field approxi- 
mation, where the homogeneous "magnetic" field B is proportional to the average order 
parameter: 

B = W((a)). (3) 

The double brackets denote both the thermal and the random-field average. 

To end the explanation of our model, a few remarks are in place. First, by replacing the 
dimerization by a discrete Ising variable, we have neglected the true dimerization profile 
associated with a soliton. As the extent of this profile is given by the correlation length £ 0) 
our " sudden- kink approximation" is valid as long as the soliton density is small compared 
to l/£o, as is the case for weak disorder. Taking into account the true dimerization profile 



results in an effective increase of the kink creation energy of the order y e£o/ a (with a the 
average carbon-carbon distance in trans-polyacetylene) and, thus, in a small reduction of 
the soliton density. This has recently been confirmed explicitely in numerical simulations 
which do account for the true profile f2l| . Second, it should be kept in mind that the 



RFIM Eq. ([I]) is an effective model, obtained by integrating out small lattice fluctuations. 
As a result, the kink creation energy /i weakly depends on the temperature ||. Third, 
above we have not specified the value of the Ising lattice constant d (which should not be 
confused with the lattice constant a of the polymer chain). It should be noted that both 
the disorder strength e and the interchain interaction energy W scale proportional to d. 
In Sec. [|we show that all physical observables are ^-independent in the d — > limit. 

The temperature versus disorder strength phase diagram of the model Eq. ([I]) contains 
two phases: the ordered phase characterized by a nonzero value of the average dimer- 
ization, ((a)) Ao, and the disordered phase, in which the long-range bond order (LRBO) 
is destroyed by thermal and disorder-induced kinks. The two phases are separated by a 
second-order transition. Figure 1 shows the phase diagram calculated for W/ (J, = 0.008, 
a typical value for trans-polyacetylene if the interchain interactions are dominated by 
interchain electron hopping [p2|| . The stars in Fig. 1 denote the phase boundary which 
we obtained by numerical simulation of the model Eq. (|I]) using an algorithm based on 
the transfer- matrix approach (cf. Ref. 0). The order parameter ((a)) was found from a 
self-consistent calculation of the mean field B and the critical curve was then obtained 
by requiring that ((a)) — ► 0. The smooth temperature dependence of the phase bound- 
ary was obtained by averaging the free energy over 10 4 random-field realizations for a 
chain with 10 3 sites. The solid curve in Fig. 1 indicates the phase boundary which was 



calculated in Ref.|18[ from an analytical expression for the average free energy of the 
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continuum version of the model Eq. (JJ). With the exception of a small temperature 



region, T < T (e) |18fl , the results of the continuum and the discrete models agree well. 

At low temperatures the phase transition results from a competition between the 
disorder and the interchain interactions. In fact, the critical disorder strength which 



separates the phases with and without LRBO at zero temperature reads |TB 



e c « Iwp . (4) 

In the disordered phase, ((a)) = 0, the density of disorder-induced neutral solitons 
(spin-flips) is to lowest order in e given by |L8] 



Us = \j^ (5) 

as is the case for a single disordered chain (W = 0) [p^j . On the other hand, for e < e c the 
order parameter ((a)) is observed to increase rapidly [|Tj| with a slope that is proportional 
to the ratio fi/W ^> 1. Thus, in the overwhelming part of the ordered phase the system 
is nearly perfectly ordered with an order parameter close to unity, ((a)) ~ 1, and the 
solitons are bound into pairs by the interchain interactions. Well within the ordered 
phase (e <C e c ) their density is exponentially suppressed. The distance between the 
soliton-antisoliton pairs is much larger than the typical pair size and the number of 
soliton-antisoliton pairs per unit length reads 



12W 2 ( Wp 



n p = - — — exp ( - 2 ——J (6) 

(for W 2 < e < 2W>/3). 

In the following we will focus on the LRBO phase and calculate the magnetic suscep- 
tibility due to the spins of the bound pairs of neutral solitons. 



3 Magnetic susceptibility in the ordered phase 

Apart from the interchain interaction discussed in the previous section, there is also an 
intrachain interaction between kinks. The latter interaction is strong only when the 
distance between kinks is of the order of their size, £q. Thus, for weak disorder, when the 
density of kinks is small, it has little effect on the statistics of the kinks. It may, however, 
be important for the magnetic properties of disordered Peierls systems, as it results in an 



antiferromagnetic exchange between the spins of neutral kinks [17]. This exchange can 
bind the spins of neighboring kinks into nonmagnetic singlets, thus reducing the magnetic 
susceptibility of the system. 

As argued in the previous section, for a nearly perfectly ordered system the typical 
distance between disorder-induced soliton-antisoliton pairs is much larger than the typical 
pair size. We may then neglect the spin exchange between kinks from different pairs. The 
Hamiltonian describing the interactions of soliton and antisoliton spins, Si and 5*2, within 
one pair reads 

H pair = j^.^-i)- gfiB H (St + S z 2 ) , (7) 
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where J is the exchange constant in the pair and H is the external magnetic field. The 
pair free energy is given by 



f{J,H) = - J - - In [l + e~ pJ (1 + 2 cosh{(3gp B H)) 
P 

and the zero-field magnetic susceptibility of the pair is 



OH 2 



e 



-pj 



= 2g^\fl - (9) 

H=0 



The coupling J decreases with the soliton-antisoliton separation R. Quite generally, 
the large- R behavior is 

J = Jo exp (-^j , (10) 

where p = ^o/d (R is measured in units of d) and Jo is of the order of the spin gap. For 
the SSH model, in which Coulomb interactions are neglected and the spin gap equals the 



charge gap, J = 4A [[L7j . 

As R is a random quantity that is imposed by the disorder realization, also J is 
random. If we know the pair size distribution, p(R), the distribution of exchange values, 
w(J), can be obtained using Eq. ([TO). We normalize the latter to the total density of 
spin-pairs 

roo 

n p = / dJw(J) . (11) 



The system's magnetic susceptibility is then given by 

X (T) = / dJw(J) X (T,J) , (12) 



with x(T, J) as in Eq. (|). 

Clearly, the temperature dependence of the magnetic susceptibility is determined by 
the pair size distribution. As we will show in detail in Sees. |] and |5|, in the LRBO phase, 
p(R) is sharply peaked at some R*, while for R^> R* 

p(R^>R*) oc exp^-a^ (13) 

with a a constant determined by the strength of the disorder and interchain interactions. 
Equation (|T3"D in a straightforward way yields a power-law exchange distribution 

w(J) oc (^)' (14) 

for J <^ J(R*). This part of w(J) dictates the behavior of the magnetic susceptibility 
at low temperature, T < J(R*). Pairs with J > J(R*) (or: R < R*) are in the 
nonmagnetic singlet state at these low temperatures. We thus find 

7 \ 1— « 

Jo 



X (T « J(R*)) oc ( ^ ) , (15) 
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which deviates from the high temperature Curie behaviour. 

We note that, in order to describe the anomalous temperature dependence of the 
magnetic susceptibility of charge transfer salts, Clark et al. [53] also introduced pairs of 
spins with a random antiferromagnetic coupling. These pairs, however, were introduced in 
a purely phenomenological way, whereas in our model they naturally emerge as disorder- 
induced soltion-antisoliton pairs with a distribution of exchange constants that follows 
from the pair size distribution. 



4 Calculation of the Pair Size Distribution 

The pair size distribution p(R) is defined as the number of soliton-antisoliton pairs of 
size R per site of the Ising chain. For a given disorder realization {h m } one finds from 
Eq. (|IJ) that the energy of a configuration (mi,m 2 ), with the soliton located between m\ 
and mi + 1 and the antisoliton between m 2 and m 2 + 1, reads: 

E[m 1 ,m 2 ] = Eq - AE[m 1 ,m 2 \ , (16) 

where Eq denotes the energy for a configuration without solitons and 

AE[m 1 ,m 2 ] = -2/i - 2 ]T (B + h m ) (17) 

m=mi+l 

is the energy change due to the creation of the soliton-antisoliton pair. As we restrict 
ourselves to isolated soliton-antisoliton pairs, it is sufficient to consider a segment of the 
chain which contains one such pair located far away from its end points. Furthermore, 
because the sequence of solitons and antisolitons along the chain is determined by fixed 
boundary conditions for the lattice dimerization, we may, without loss of generality, 
assume that m 2 > m-y. Then, the pair size R (in units of the Ising lattice constant d) is 
given by 

R = m 2 — mi . (18) 
The soliton-antisoliton pair configuration (mi,m 2 ) is only energetically favorable if 

AE[m 1 ,m 2 \ > . (19) 

This is, however, not sufficient to calculate the pair size distribution p(R), as we also 
have to impose the condition that this pair configuration has lower energy than any 
other pair in the considered chain segment. Thus, simultaneously, the energy of the pair 
configuration has to satisfy the inequalities 

AE[m 1} m 2 ] > AE[m Xl m' 2 ] , (20) 

for all other possible pair configurations (m^m^). Therefore, the desired pair size distri- 
bution takes the form: 

p(R) = /e(AE[m 1 ,m 2 \) U 6 ( AE[m u m 2 ] - AE[m' x , m 2 ] )\ , (21) 
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where Q(x) is the step function: 

^ = {1 t:-<l ^ 

and the brackets, (...), denote the Gaussian average over disorder realizations {h m }. The 
definition Eq. ([H]) ensures that, in accordance with Eq. (|llD, p(R)/d is normalized to 
the density of soliton-antisoliton pairs: 



1 

d Jo 



n p = - I dRp(R) . (23) 



It is easy to see that p(R) factorizes into two independent parts, p out and p in , that 
account for the soliton-antisoliton pairs with a size that is, respectively, larger and smaller 
than R: 

p(R) = p out Pi n , (24) 

with 

Pout = ( ®(z + s mi ) Q(2z + s mi + *mj — 1 

) Q(3z + s mi + 

Q{z + s m2+ i) Q(2z + s m2+ i + s m2+ 2) Q(3z + s m2+ i + s m2+ 2 + s m2+ 3) ... ) , (25) 



and 

/ m2 \ 

p in = ( e(-j -Rz- Y, s m) n L n R ) . (26) 

m=mi+l 

Here, we have defined the dimensionless variables s m = h m /\/e, z = B/y/e, I = n/y/e, 
while 

n L = Q(-z - s mi+ i) 6(-2^ - s mi+ i - s mi+2 ) ... ®{-Rz - s mi+1 - ... - s m2 ) (27) 
and 

n fl = ©(-2 - s m2 ) Q(-2z - s rn2 - s m2 _i) ... Q(-Rz - s m2 - ... - s mi+1 ) . (28) 

Note that p out itself also consists of two independent factors: the first factor excludes 
the pairs with the soliton located to the left of mi, while the second one excludes anti- 
soliton positions larger than rri2 + 1. Both these factors can be written in the form: 



y{z) = n 



m=l 

where 



/+oo , m 

ds m f(s m ) 9( + S k 

■°° k=l 



(29) 



1 „2\ 



V Z7T 



is the Gaussian weight. As a result, for the outer factor we obtain: 



Pout = [Y(z)f . (31) 
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The function Y(z) will be calculated later in this section. 

The calculation of the inner factor is complicated by the presence of the extra 0- 
function in Eq. (p6[), which precludes the factorization of pi n in two independent averages. 
However, considerable simplification is possible, because we focus on the bond-ordered 
phase where ((er)) — 1. Then the density of disorder-induced soliton-antisoliton pairs is 
small, and the main suppression factor in p(R) is the probability of the disorder realization 
necessary to create a pair PR In other words, the most important contribution to 
Pin [and p{R)\ comes from averaging the first G-function in Eq. (f26|): 



11(2 



e(-/- £ (z + sm))) = leMg(R)) « , (32) 



m=mi+l 

where 

a(R) = —, 

Here, the asymptotic expression for the complementary error function erfc(g(i?)) was 



g(R) = ^£ • (33) 



used because the minimal value for its argument is easily shown to be g m i n = y3e c /e, so 
that for e < e c /2 the relative error becomes already less than several percent. 

The interpretation of this result is that the optimal disorder fluctuation (i.e., the 
disorder realization with the largest weight) that can induce a soliton-antisoliton pair 
of size R has a constant value —Hr in the interval of length R and is zero outside the 
interval. The amplitude h R is determined by the energy balance [see Eqs. ([P7D and (|19|) 1: 

h R R = fj, + WR . (34) 

The weight of the optimal fluctuation, w = exp (^—RhR 2 /(2e)^j, is precisely the exponen- 
tial factor appearing in Eq. (|32J). At 

= £ = (35) 

the weight reaches its maximal value, exp (—2Wfi/e). For e e c , the maximal weight 
is small [as was also found in Eqs. ( |3~2l) and (|3~3D 1 and the soliton-antisoliton pairs are 
suppressed. In that case, all disorder realizations that contribute significantly to p(R) 
are close to the optimal fluctuation. 

Bearing this in mind, we now calculate the inner factor Eq. First, we can rewrite 
Eq. ( p6|) in the form: 



Pi,. 



poo r+ioorj\ " 1 2 f+oo 

/ dS / — e- xs J] / ds m e- Xs -f(s m )U L U R , (36) 

JzR+I J -too Z7TZ m=m 1 +l 00 

where the integration over A ensures that 

m=m,2 

S = - £ s m (37) 

m=mi+l 



and the limits of the integration over S follow from the first G-function in Eq. p 
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Since the typical pair size R* ^> 1 [see Eq. ([35])], we can use the canonical formalism, 
in which Eq. (|37|) for the sum of R random variables is satisfied only in average. We do 
this by "shifting" the argument of the random-field distribution on each site by A: 

/(s ) _> f {s + A ) = e -^ 2 - Xs f(s) , (38) 
so that the average value now becomes s = —A and Eq. ([JBp reads: 

Pin = dS e- xs+ * RX J] / ds m f(s m + \)U L U R . (39) 

JzR+I J-ioo Z7U m=mi+l 

The integral over A comes from the small vicinity (~ 1/Vr,) of Ao = S/R, where the 
exponential in Eq. (|59"|) has its maximum. Saddle-point integration over A then gives: 



Pi 



dS 



R+i V2nR 



e 2 « 



™2 „ + 00 / 

II / ds m f [S m + — ) r^ITR 

= mi+ l J ~ 00 V K/ 



(40) 



Next we note that if the condition imposed by the first G-function in Eq. (|2BJ) is 
satisfied, the arguments of the last 0-functions in and 11^ also almost certainly are 
positive. In other words, because the relevant disorder realizations are close to the optimal 
fluctuations, only the first few G-functions in 11^ and Il# are really restrictive. This 
implies that the disorder averages of 11^ and Ur in Eq. ( |4"0"D are decoupled. Furthermore, 
it is easily seen from Eq. (g|) that then (II L ) = (Tl R ) = F(| - z), so that Eq. (g^) 
becomes: 



Pi 



dS 



-.e 2 « 



izR+i \ / 2tcR 

The integral over S comes from the vicinity of the lower limit, S = zR + I [cf. Eq. (|3 
for the optimal fluctuation]. The result of the integration is: 



(41) 



R exp 



Pin \ „ 




(I+Rzf 



2R 



2tt (/ + Rz) 



(42) 



where for S in the argument of Y we took its value at the lower limit of the integration. 
From Eqs. (|24]), (|3l|), and (^2|) we finally obtain for the pair size distribution: 



p(R) 



exp(-g(R)< 



^g(R) 



-4 



[Y(z)f 



(43) 



where the function g(R) is defined by Eq. (|33|). 

What is left now, is the calculation of the function Y(v ). To this end we introduce the 
function Y(s\v), satisfying the integral equation: 



Y(s\v) 



ds' f{s + v-s')Y{s'\v) . 



Comparing the iterative solution of this equation to Eq. (|29|), one finds: 

Y(v) = Y(0\v) . 



(44) 



(45) 
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The integral equation (|4l ) can be easily solved numerically. The result is shown as stars 
in Fig. 2. The solid line represents the best fit to these points by a function of the form: 



Y(v) = tanh(cf) . (46) 

The fit yields c ~ 1.14. For small v, the best linear fit (Y(v) = c'v) yields d ~ V2, with 
a precision of several percent. 

We conclude this section by a brief analysis of the pair size distribution Eq. (|43|). First, 
integrating p(R)/d over R gives the density of soliton-antisoliton pairs [cf. Eq. (pffi) ]. For 
e <C e c , the exponential factor in Eq. Q4"3| ) has a sharp peak at R* given by Eq. (|35|). 
Using B ~ W, the saddle-point integration around the peak gives: 



n 



v d2W 2 



Y 



(5 



exp -2— . (47) 



Furthermore, we note that knowledge of the pair size distribution Eq. fl43|) allows us to 
derive the long- and short-range bond order parameter of the RFIM Eq. ([[]). For a dilute 
gas of soliton-antisoliton pairs it is sufficient to consider a chain segment of iV + 1 sites 
containing a single pair of size R: 

a(m) = 1 — 2 0(m — rrii) 0(m 2 — m) , (48) 

where mi and ui2 = mi + R denote the positions of, respectively, the soliton and the 
antisoliton. Replacing summations by integrations, the LRBO parameter averaged over 
all possible pair sizes is easily calculated: 

((a)) = / dRp(R) / n dma{m) = 1 - 2 / dRp(R) R . (49) 

J J 2~ 

For a nearly perfectly ordered system, ((cr)) ~ 1, we thus find from Eq. ( p9|) the con- 
dition that the typical pair size R* is much smaller than the typical number of sites 
l/(dn p ) between soliton-antisoliton pairs. Similarly, we calculate the correlation function 
((<t(0)<t(Z))), which yields the sum of the square of the LRBO parameter Eq. (f49"D, 



roo 

((a)) 2 ~ 1 - 4 / dRp(R) R , (50) 
Jo 

and the connected correlator 

<M0)<7(Z)» C = 4 / dRp(R+\l\)R. (51) 
Jo 

The scale for the decay of short-range correlations is obviously set by the typical pair size 
R*, as ((a(0)a(l))) c -> for \l\ > R*. 



5 The Exchange Distribution in the Continuum Limit 

In the previous sections we described disordered Peierls systems using the effective RFIM 
Eq. (jl]). The values of the interchain interaction W and the disorder strength e in this 
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model are proportional to the length d which we choose for the unit cell of the Ising 
chain and which plays the role of a short-distance cut off. On the other hand physical ob- 
servables, such as the density of soliton-antisoliton pairs and the magnetic susceptibility, 
should not depend on d. Thus, before comparing our results to the experimental data on 
irans-polyacetylene, we show that d drops out from the expressions for the observables 
in the d — > limit. 

To this end, we introduce as physically meaningful quantities the disorder strength e 
and the interchain interaction W per unit length: 

e = ed, (52) 

and 

W = Wd. (53) 

Furthermore, from now on we will work with the physical pair size r = Rd. In terms 
of these new variables, the arguments of both y-functions in Eq. (JE|) for the pair size 
distribution are cx \fd. Therefore, for d — ► 0, the arguments are small and we can use 
Y(v) ~ V2v [see below Eq. giving the pair size distribution 

p{r) = 4f w»vpm (54) 

e z r z y 47r gir) 

where 

/i + Wr 

g V> = — 7k=^~ ■ ( 55 ) 
V2er 

Similarly, from Eq. fl47|) , the total density of neutral soliton-antisoliton pairs in the limit 
d — > is found to be: _ 

n p « - — exp [-2^ ] , (56) 



which coincides with Eq. (||) obtained in Ref. []I8] within the continuum approximation 
for the RFIM Eq. (0). 

Furthermore, in terms of continuum variables, the exchange coupling Eq. fllHD reads: 

J(r) = Jo exp(-f ) . (57) 
Thus, in the continuum limit the distribution of exchange constants becomes: 

w(J) = / drp(r)5(J-J(r)) = 4 - - p T » 58 

Jo J e 2 [\n(J/J )\ 2 V47T p( J) 

with <7(J) = g(r = e O ln(J /J)) [cf. Eq. ©]. 

In Fig. 3 we plot the distribution w(J) for four different parameter sets (e, W, Jo) 
chosen such that the density of soliton-antisoliton pairs is fixed at n p = 1/6000 a -1 (with 
a the average carbon-carbon distance in trans-polyacetylene). Our choice of parameters 
is summarized in Table I and will become clear in Sec. [| Depending on the parameters, 
one observes two qualitatively very different behaviors: w(J) either has a pronounced 
peak at J* ~ J(r*) (with r* = R* d and R* as in Eq. fl35D ) and tends to zero for J — > 0, 
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or w(J) diverges for small J. The distinction between these two behaviors is dominated 
by only one parameter combination: 



In fact, Eq. (EH) for J < J* yields 



W 2 

a = — Co. (59) 



W (J) cc (y) 1 ° , (60) 

which shows that the relative strength of the interchain interactions and the disorder 
determines whether w(J) diverges (a < 1) or approaches zero (a > 1) for J — > 0. 

The behavior of Eq. ( ptj) agrees with Eq. flTj) in Sec. [3] and can indeed be traced back 



to the fact that for large r the pair size distribution is exponential: 

( ^ *\ ( w " r \ ( J °Y a (au 

p[r »r ) oc exp — - =1—1 . (blj 



2e J \J 

This exponential dependence can be understood as follows. For r ^> r*, the energy 
of the string between soliton and antisoliton exceeds the kink creation energy: Wr = 
WR ^> fi. Thus, the amplitude of the optimal fluctuation Eq. (|34| ) is HrR ~ WR. 
The Gaussian weight exp (— Rh 2 R /(2ej) of this fluctuation is the exponential in Eq. (|6lD . 
Similar arguments were used to explain the power-law dependence of the density of states 
in the Fluctuating Gap Model of disordered systems pB]] . The power-law dependence of 
w(J) at small J gives rise to a characteristic low-temperature behavior of the magnetic 
susceptibility, as we will see in the next section. 



6 Magnetic susceptibility of £ran5-polyacetylene 

In this section we consider irans-polyacetylene as a disordered Peierls system and calculate 
its magnetic susceptibility due to disorder-induced soliton-antisoliton pairs as a function 
of temperature. The temperature dependence of the magnetic susceptibility is determined 
by the distribution of exchange constants Eq. (|58|) . For temperatures T much larger than 
the typical singlet-triplet energy splitting J*, almost all spin pairs are thermally excited. 
Thus, we have, essentially, 2n p free spins, which give rise to a Curie susceptibility. Indeed, 
Eqs. (|) and (|IJ) yield: 

1 r°° 1 7i 

X (T » J*) ~ -<7 2 /i|/3 J Q dJw(J) = -g 2 ^ ^ , (62) 

where the density of soliton-antisoliton pairs n p is given by Eq. (p6|). In the opposite limit, 
T ^ J*, however, most of the spin pairs are in the singlet state and do not contribute 
to the magnetic susceptibility. Under these conditions Eqs. dl2|) and fl58D yield: 

X (T « J*) = C{T) , (63) 



13 



with a as in Eq. (|59"D and logarithmic temperature corrections given by 

cm ~ i^i (^k\ h m 2 r(i + a )Li a (-|) 

1 J ~ 9J « ^ vr ; e£ 2 (ln(J /T))V2 
(Li„(;z) = Y^k=\ zk is the polylogarithm function). The dominant factor in Eq. (|6 



is (Jo/T) 1 a , which basically gives the density of spin pairs with singlet-triplet splitting 
~ T [cf. Eq. ©]. 

We thus find that the low-temperature behavior of the magnetic susceptibility differs 
from the Curie law and is dictated by the relative strength a of interchain interactions and 
disorder. For a < 1, the susceptibility diverges as T — > 0, while for a > 1 it approaches 
zero. 

A low-temperature (T < 30 K) deviation from Curie behavior has indeed been ob- 
served by Foot et al. in ESR experiments on Durham trcms-polyacetylene [[5J. These 
authors already suggested pairing of spins as possible reason for this behavior. To see 
whether our model of spins associated with disorder-induced soliton-antisoliton pairs of- 
fers a microscopic explanation, we compared our result for x(T) [numerically calculated 
from Eqs. (|T2|) and fl58)] to the experimental data. In our fit procedure there are, in 
principle, three free parameters: W, e, and Jo (f° r the correlation length we take the 
SSH value £ = 7a). We require, however, that our parameters are also consistent with 
the total density of spin pairs, which is reported to be approximately n p = 1/6000 a -1 
@? H> B 01- This requirement imposes a relation between W and e, reducing the number 
of free parameters to two. 

In practice, we chose various values for a [Eq. Q59p]. For each a value, W and e 
are uniquely determined by n p , and J is left as free parameter to fit the temperature 
dependence of the magnetic susceptibility. This procedure yielded the fits shown in Fig. 4, 
with parameter sets given in Table I. For convenience, we will refer to each parameter 
set by its a value. We note that our values for J are much smaller than the value 
Jo = 4Ao ~ 10 4 K, which one would expect in the absence of Coulomb interactions. It is 
not known, however, how Coulomb interactions alter Jo and (possibly) the exponent in 
Eq. (0). 

The important point is now that, while all four parameter sets give rise to reasonable 
fits of the experimental data, they predict totally different behaviors for T < 5 K, where 
experiments have not been performed. This is shown in Fig. 5, where we extend the 
four theoretical fits to 1 K. The qualitative differences in the low-temperature behavior, 
dictated by the value of a, are clearly visible below 5 K. This suggests that extending the 
experiments to lower temperatures may yield more information on the relative strength 
of interchain interactions and disorder in trans-polyacetylene. 

We conclude this section by noting that, within the context of our model, it is possible 
to determine the strength of the interchain interactions and the disorder independent of 
the fitting parameter Jo. For this purpose, the density n p of soliton-antisoliton pairs is to 
be obtained experimentally from the Curie tail of the magnetic susceptibility, while for 
the same sample a is to be determined from the asymptotic zero-temperature behavior of 
the magnetic susceptibility. Then, using Eqs. (|56|) and (59|), the strength of the interchain 
interactions, 



W ss 4a — 



(65) 
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and the disorder strength, 



In 



4a 



(66) 



can be calculated function of n p and a. 



For typical values of the density n p the 
logarithmic factor depends only weakly on a ~ 0(1) and can be approximated by a 

7a, 



numerical constant. If n p = 1/6000 a 



and choosing again the SSH-parameter £ 
0.02 a fi 2 / a. It is important to realize that Eqs. (65) 



one obtains W ~ 0.07au/aande 
and (B^) do not depend on the maximal exchange constant Jo which may be used as a 
fitting parameter for temperatures around the typical exchange J*. 



7 Concluding remarks 

To summarize, we have calculated the magnetic susceptibility of quasi-one-dimensional 
Peierls systems with a doubly degenerate ground state. We have related the temperature- 
dependent part of the susceptibility to the presence of neutral solitons and antisolitons 
with spin | induced by disorder in the electron hopping amplitudes along the chain. We 
have assumed the interchain interactions to be sufficiently strong to bind the disorder- 
induced solitons and antisolitons into pairs and thus establish long range bond order in 
the system. Using a mapping on the random field Ising model, we have calculated the 
distribution of the soliton-antisoliton pair size. This allowed us to obtain the distribution 
of exchange constants describing the interaction between the spins of the soliton and the 
antisoliton within one pair. Both distributions strongly depend on the relative strength 
a of the disorder and the interchain interactions. As a result, the magnetic susceptibility 
deviates from the Curie law: below T = J*, where J* is the most probable value of the 
exchange constant, the magnetic susceptibility behaves as (1/T) 1_Q . 

Our results explain the deviation from Curie behavior observed in Durham trans- 
polyacetylene though from the experimental data it is difficult to find unambigously 
the values of a and J* for this conjugated polymer. It is, therefore, important to ex- 
tend the measurements to lower temperatures, where the temperature dependence of the 
magnetic susceptibility is extremely sensitive to the choice of parameters. 

Our theory is only applicable when the low-temperature behavior of the susceptibility 
is an intrinsic property of the material and is not governed by spins of impurities. The 
latter situation may, in fact, be realized in Shirakawa trans-polyacetylene, which shows 
Curie behavior down to T = 1.5 K [E7J. Furthermore, we assumed the existence of long 



range order in the system. Whether this is the case in trans-polyacetylene is an open 
question. It would, therefore, be interesting to extend our studies to the case without 
long-range order. At the same time, however, it should be noted that in substituted 
polyacetylenes, the degeneracy of the two dimerized configurations may be lifted. This 
leads to an extra (intrachain) source of soliton-antisoliton confinement p^] and favors 
long-range bond order. Our theory may be applied to these substituted polymers by 
simply adding to the interchain interaction per bond (Wa), the energy difference per bond 
between the two dimerized configurations. As this energy difference may be controlled by 
varying the substitutions, this opens interesting possibilities to study disorder-induced 
solitons in more detail. 
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We finally mention that in the ordered phase, the disorder-induced soliton-antisoliton 
pairs show up in the x-ray spectrum as a broad incoherent peak associated with each 
sharp elastic peak arising from the bond length alternation. Our result for the pair size 
distribution allows one to calculate the shape of this incoherent peak: it simply is the 
Fourier transform squared of the connected correlator Eq.(|5"T|). Thus, one immediately 
finds that the peak width is ~ 1/r* ~ W/fi. It should be kept in mind, however, that 
this calculation does not account for other broadening mechanisms, e.g., those due to the 
complicated morphology of polyacetylene samples. 
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Tables 



W {n/a) 


e (fx 2 /a) 


Jo (K) 


a 


0.038 


0.010 


400 


0.5 


0.055 


0.014 


130 


0.75 


0.070 


0.017 


72 


1.0 


0.100 


0.023 


39 


1.5 



Table I. The four sets of parameters used in the numerical calculations for a fixed density 
of spin pairs, n p = 1/6000 a -1 . The value for a as defined in Eq. fl5"9|) is obtained using 
the SSH-parameter £ = 7a for the correlation length. 
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Figures 



vu 3 




Fig. 1. The phase diagram of the RFIM Eq. ([!]) captures the essential physics of weakly 
disordered Peierls systems and is shown as a function of the disorder strength e//j, 2 and the 
temperature T j \i. The long range bond order (LRBO) phase corresponds to an average 
dimerization ((a)) Aq ^ 0. The numerical calculation of the critical curve (stars) agrees 
well with the analytical result (solid curve) obtained in Ref. fl8| . The dashed curve 
indicates the breakdown of the continuum approximation in the analytical calculation 
below T = T (e) W 
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V 



Fig. 2. The function Y(v) obtained from a numerical solution (stars) of the integral 
equation d44[) . The best fit of this solution by a function of the form Y(v) = tanh(cv) 
yields c ~ 1.14 (solid line). 
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Fig. 3. The distribution w(J) of exchange constants as a function of J/J for the four 
parameter sets given in Table I. The curves correspond to a — 1.5 (dots), 1.0 (dash-dot), 
0.75 (dashes), and 0.5 (solid). For a > 1.0 the distribution has a pronounced peak at 
some J = J* and tends to zero for J — > 0. In contrast, for a < 1.0, w(J) diverges when 
J ^ 0. 
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10 20 10 20 
100/T (1/K) 100/T (1/K) 



Fig. 4. Fits of our theory (solid curve) to the experimental data (dots) for the magnetic 
susceptibility of Durham trans-polyacetylene obtained in Ref. 0. The four parameter 
sets given in Table I were used to fit the same experimental data points (see text for 
details). The deviation from Curie behavior (straight line) below T = 30 K is clearly 
seen and reasonably reproduced by each fit down to T = 5 K, below which experimental 
data are not available. 
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a = 1.5 



20 40 60 80 100 
100/T (1/K) 

Fig. 5. As Fig. 4, but now the four theoretical curves are shown down to T = 1 K. 
It is clearly observed that different values for a (the relative strength of disorder and 
interchain interactions) lead to qualitatively different low-T (< 5 K) behavior of the 
magnetic susceptibility. 
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